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Abstract. This paper presents the first I-band photomet- 
ric catalog of the brightest galaxies extracted from the 
Deep Near Infrared Survey of the Southern Sky (DENIS) 
An automatic galaxy recognition program has been de- 

' veloped to build this provisional catalog. The method is 
based on a discriminating analysis. The most discriminant 
parameter to separate galaxies from stars is proved to be 
the peak intensity of an object divided by its array. Its 
efficiency is better than 99%. The nominal accuracy for 
galaxy coordinates calculated with the Guide Star Catalog 

' is about 6 arcseconds. The cross-identification with galax- 
ies available in the Lyon-Meudon Extragalactic DAtabase 
(LEDA) allows a calibraton of the I-band photometry with 
the sample of Mathewson et Al. Thus, the catalog contains 
total I-band magnitude, isophotal diameter, axis ratio, po- 
sition angle and a rough estimate of the morphological 
type code for 20260 galaxies. The internal completeness 
of this catalog reaches magnitude lu m — 14.5, with a pho- 
tometric accuracy of ~ 0.18m. 25% of the Southern sky 
has been processed in this study. 

This quick look analysis allows us to start a radio and 
spectrographic follow-up long before the end of the survey. 

Send offprint requests to: G. Paturel 

* Based on observations collected at the European Southern 
Observatory, La Silla, Chile. 
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1. Introduction 

In Paturel et al. (1996), we presented our program of col- 
lecting the main astrophysical parameters for the principal 
galaxies. The first target was limited to adding informa- 
tion on galaxies already known in the LEDA database. 

The work is now more ambitious because we are aim- 
ing at detecting new galaxies from the Deep Near In- 
frared Survey of the Southern Sky (hereafter DENIS). DE- 
NIS is a program to survey the entire southern sky in 
three wavelength bands (Gunn-i: 0.80//to, J: 1.25/im and 
Ks: 2.15/xm) with limiting magnitudes of 18.5, 16.5 and 
14.0 respectively. The observations are performed with 
the ESO lm-telescope at La Silla (Chile) with a dedi- 
cated camera. The survey observations with three chan- 
nels started in routine mode in December 1995. A detailed 
description of DENIS is given in Epchtein (1998) and in 
Garzon et al. (1997). 

The systematic detection, extraction and cataloging 
of DENIS extragalactic sources are of significant inter- 
est for studies requiring large and homogeneous samples 
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such as the kinematics of the local universe, the distance 
scale, cosmology etc... The I-band is the most suitable 
both for the detection of extended objects and for the 
star/galaxy separation, except in the galactic plane. Thus, 
in the present study only the I-band images are considered. ■ 9 ' 

The transfer of DENIS images from Paris to Lyon is * 
explained in section 2 and the extraction of objects from 
these images in section 3. In sections 4 to 5 we describe as- 
trometry and automatic galaxy recognition and analysis. 
Then, in sections 6 to 8, we explain how galaxies are cross- 
identified with LEDA galaxies leading to the comparison 
of astrophysical parameters with those from Mathewson 

et al. and from LEDA. Finally, in section 9 we describe ' ■ *■ ■ 

the provisional I-band DENIS catalog. 

It is important to note that a deeper catalog will be 
made at the end of the survey. Hence, the present cat- 
alog is a preliminary catalog of bright galaxiesdetected 
by DENIS, and used to start a radio and spectrographic 
follow-up long before the end of the survey. The present ^ t 
measurements cover one year of observation. 



* 



2. Obtaining I-band CCD frames from DENIS 

2.1. Characteristics of the DENIS survey 

The I-band of the survey is the Gunn-i band at 0.80/im. 
The CCD camera is a Tektronix 1024x1024 pixel array 
cooled down to liquid nitrogen temperature. Each frame 
(768x768 pixels) represents 12 square arcminutes with a 
pixel size of 1 arcsecond. The integration time is 9 seconds. 
The read-out noise is about 7e~ . The observing strategy 
consists in scanning at a constant right ascension on strips 
of 30 degrees in declination (180 frames per strip) taken 
in three zones "Equatorial" from 5 — +2deg to — 28deg, 
"Intermediate" from S = — 28deg to — 58deg and "Polar" 
from S = —58 cleg to — 88deg. The overlap between adja- 
cent frames is 1 arcminute on each side (i.e. 2 arcminute 
in each direction). This strategy aims at covering a wide 
range of airmasses. Each strip starts and ends with pho- 
tometric and astrometric calibrations. At the end of some 
nights a flat fielding is performed directly on the sky dur- 
ing sunrise. The data is archived on DAT cartridges which 
are send each week to the Paris Data Analysis Center 
(PDAC) at the Institut d'Astrophysique de Paris for pro- 
cessing. 

2.2. Pipeline Lyon-PDAC 

A systematic automatic processing of DENIS I-images be- 
gan in Lyon in February 1996. Flat-fielding and de-biasing 
are made at PDAC on each genuine frame. For the Lyons 
processing each I-image is reduced by a factor 4 in size 
by rebinning pixels 2x2. Our effective pixel size is thus 2 
arcseconds. An example of an I- frame is given in Figure [l| 
The histogram of pixel intensities is used to calcu- 
late the sky background intensity fb g . The sky-background 
level is taken at the maximum intensity of the histogram 



Fig. 1. A typical image from the DENIS Survey after the 2x2 
rebinning. The frame is 12' x 12'. North is on the bottom side, 
East is on the left side. Among the four galaxies clearly visible 
on the top, one (the third from the left) is a new one. This 
shows that relatively bright galaxies can be discovered. 

(i.e. the mode) and the standard deviation a is calculated 
by symmetrizing the low intensity part of the histogram 
with respect to the mode. The sources are concentrated 
in the high intensity part of the histogram. 

A threshold is then applied (the threshold level is 
fbg + 2. Oct). Using this procedure (averaging and thresh- 
olding) allows a compression factor of 20 to 30, depending 
on the image contents. All images of a given I-strip are 
thus compressed, tar'd and automatically transferred to 
Lyon via ftp. A full strip is stored in 10 to 13 Mbytes. 
The galaxy extraction is made at Lyon using the program 
described in the following section. 

3. Extraction of astronomical objects 

All sources (stars, galaxies, defects etc..) are extracted 
using the same algorithm as described in Paturel et al. 
(1996, section 3.1), except that no attempt was made to 
share interacting objects which are simply flagged after 
visual inspection (section 5). The reason is that we are 
interested first of all in well defined objects. At the end 
of this stage, we obtain for each frame a collection of ma- 
trixes (see an example in Figure^). Matrixes smaller than 
17 pixels are rejected. They correspond, to the mean, of 
objects of 6 arcseconds in diameter. 

Astrophysical parameters arc extracted for each ma- 
trix according to Paturel et al. (1996, section 3.2). These 
parameters are the following: 
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Fig. 2. Example of a matrix for a small 
2// x 2//. The edge is outside the matrix. 



laxy. One pixel is 



10 20 30 

Number of GSC stars 



— Weighted mean position of the center x m and y m (in 
pixels), where the weight is the pixel intensity /y 

— major axis logD (D in O.larcmin.) at faint isophotal 
level (external diameter). 

— axis ratio logR (log of major to minor axis) 

— position angle of the major axis [3 (in degrees, counted 
from North towards East) 

— magnitude I = — 2.5Zog(^™ =1 — fb g ) + est (in arbi- 
trary units), n is the number of pixels with intensity 
fij larger than the threshold. 

We have now to perform astrometry (conversion of pix- 
els positions to right ascension and declination) and then 
recognition of "galaxies" , "stars" , and "unknown objects" . 

4. Astrometry 

The center of the frame is taken from the header of 
the FITS file. From the Guide Star Catalog (GSC) and 
from the LEDA database we extract all objects (stars 
or galaxies) known in the corresponding square. A cross- 
identification between matrixes and stars is made exactly 
as described in Paturel et al. (1996, sections 3.3 and 3.4). 
Galaxies are also used but only when they have accurate 
coordinates (i.e. typically better than 10 arcsecond). A 6th 
order polynomial fit converts (x,y) positions on the frames 
to Right Ascension and Declination. The number of GSC 
stars varies from one frame to another. A histogram of 
number of GSC stars per frame is given in Figure | If this 
number is smaller than 7 or if the standard deviation of the 
polynomial fit is greater than 4 arcseconds, the solution 
is rejected and we adopt the 'header' solution calculated 
from the coordinates of the center and the pixel size as 
given in the header. If the GSC-solution seems acceptable 
but differs from the header-solution by more than 30 arc- 
seconds, the header solution is preferred and coordinates 
are flagged to recall that they may be inaccurate. 

In Figure f|, we show the differences between the GSC- 
solution and the FITS header-solution. Most of them are 
in good agreement within 15 arcseconds. Note that more 



Fig. 3. Histogram of the number of GSC stars by frame. 



recent measurements have been astrometrically calibrated 
to better than 1 arcsecond by cross-identifying with the 
PMM database. 




10 20 30 

Shift of coordinates (") 

Fig. 4. Histogram of the differences between the two astro- 
metric solutions (GSC- and direct-solution) . 



Galaxy coordinates will be compared directly with co- 
ordinates of LEDA galaxies. 



5. Automatic star/galaxy separation 

5.1. Test sample 

A sample of 1146 objects has been visually classified into 
three classes: stars, galaxies and unknown objects. The 
distribution in each class is the following: 
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Class 


Number of objects 


percentage 


Stars 


500 


43.6% 


Galaxies 


523 


45.6% 


Unknown 


123 


10.7% 


Total 


1146 


100% 



This sample will be used as a test sample (or training 
sample) in a discriminant analysis method (DA) for galaxy 
recognition. 

5.2. Discriminant analysis 

The DA is a common method for automatic recogni- 
tion. The test sample being shared in n c classes {Ik} 
(k = l,n c ), the purpose of DA is to find the principal 
factorial axis on which each class is as concentrated as 
possible and as distinct as possible from the others. This 
is achieved by maximizing the inertia between classes and 
by minimizing the inertia within each class. Inertia is cal- 
culated from the set of parameters attached to each object. 
We will note pij the j-th parameter of object i. 

The mathematical result (see Diday et al., 1982) is 
that the factorial axes are the eigenvectors of the matrix 
T~ X .B, where T is the total covariance matrix and B is 
the inter-class covariance matrix. 

Note that the matrix T" 1 .B is not symmetrical and 
that the total covariance matrix is the sum of the intra- 
class covariance matrix W (covariance Within class) and 
of inter-class covariance matrix B (covariance Between 
class). This is called the Huyghens decomposition. 

T = W + B (1) 

The elements of B are: 

" c N 

^■'=E^(p* -p>) ( 2 ) 

fe=l 

where Nk is the number of objects in the class Ik (k=l to 
n c ), where the mean parameter j for the whole sample is: 

i 

(N is the number of objects of the whole sample), and 
where the mean parameter j within the class {Ik} is: 

The elements of the total covariance matrix T are: 

k ieh 

Now, we have to choose the set of parameters attached 
to each object. 



5.3. Choice of discriminant parameters 

Any discriminant method requires a good choice of dis- 
criminant parameters which are used for the definition of 
the metric. These parameters are not necessarily indepen- 
dent but they must cover all features which seem relevant 
for a reliable discrimination of astronomical objects. For 
galaxy recognition we tested 7 parameters. 

1. Peak intensity per area unit, this is Peak intensity di- 
vided by the surface of the considered object. 

2. Mean surface brightness, total flux divided by area 

3. Peak intensity, 

4. Axis ratio, ratio of the major to the minor axis 

5. Relative area, ratio of number of pixels of the object 
and of the matrix. 

6. Elongation of the matrix 

7. Presence of diffraction cross 

The DA method is applied on half the sample (i.e. 573 
objects) and tested on the other half using only one pa- 
rameter at a time (in this case the factorial axis is defined 
by the parameter itself). The percentage of good results 
is given below for each one, individually. 



Parameter 


stars 


galaxies 


mean 


Peak over area 


100.0% 


98.9% 


99.4% 


Mean SB 


98.7% 


99.6% 


99.2% 


Peak intensity 


96.2% 


98.5% 


97.4% 


Axis ratio 


83.7% 


54.8% 


69.2% 


Relative surface 


77.0% 


58.5% 


67.7% 


Elongation of matrix 


88.3% 


40.8% 


64.5% 


Diffraction Cross 


68.2% 


57.0% 


62.6% 



The conclusion of this test is that the most relevant in- 
formation about the nature of an object is contained in the 
pixel intensity, not in the shape of the object. Stars have 
a very high central intensity, galaxies do not. Moreover, 
stars are concentrated, galaxies are not. This explains why 
"Mean SB" and "Peak over area" give such an impressive 
recognition rate. Finally, only the first four parameters 
have been used. The axis ratio is kept because it becomes 
relevant for faint objects despite that its rate is relatively 
low. 

5.4- Result 

The DA method is applied with the four parameters de- 
scribed above and three classes "Galaxies", "Stars" and 
"unknown objects" . Using the test sample, each object is 
projected onto the first factorial axis. Figure || shows the 
projection onto the first factorial axis of "Galaxies" and 
"Stars" classes. Similar plots exists for "Stars" and "un- 
known objects" classes and for "Galaxies" and "unknown 
objects" classes. All "unknown objects" have been elimi- 
nated in the next part of this study. 

One can see that there is an overlapping region where 
"Galaxies" and "Stars" are mixed. The limits of this zone 
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can be tuned in such a way that one can accept a given 
percentage of misclassification. We choose 0% chance of 
classifying a star as a galaxy and 5% chance of classifying 
a galaxy as a star. Indeed, it is important to avoid the 
contamination of the catalog by stars while it is not as 
important to miss a galaxy (which is uncertain anyway). 
These limits are drawn on Figure || where it is visible that 
no star enter the galaxy-domain, while 5% of galaxies en- 
ter the star-domain. Objects between these two limits will 
be classified as undefined. 
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Galaxies 



Stars 



-0.2 



0.2 
First factorial axis 
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Fig. 5. Definition of acceptation zones along the first factorial 
axis. The left-hand zone defines "Galaxies", the righthand one 
defines "Stars", and the intermediate one defines "undefined 
objects" . 



5.5. Visual control 

The final step of this treatment consists in checking vi- 
sually all frames recognized as galaxies. This tedious part 
allows us to reject artefacts (1148 rejections after the in- 
spection of 54073 images) like those produced by star halos 
truncated by the edge of the frame. Such truncated halos 
look like elongated, low-surface brightness object, easily 
accepted as galaxies. 

As a result, a code is given to describe three features: 

— "multiple" , if several objects are present in the matrix 

— "truncated" , if the galaxy is truncated by the edge of 
the array. 

— "peculiar" , if the galaxy looks strange for any reason 

So, each galaxy of the catalog has been inspected visually. 
This will prevent us from gross misidentification. Now, the 
galaxies have to be cross-identified with known galaxies. 

6. Cross-identification 

To make a correct cross-identification we need coordinates, 
calibrated magnitudes and diameters, axis ratios and posi- 
tion angles. We will use the first calibration obtained from 



the preliminary cross-identification made for astrometric 
purposes (section 4). This calibration will be refined in 
next section. 

Because of frame overlap along the strip, many objects 
are measured twice (or even three or four times with ad- 
jacent strips). A first cross- identification is done for these 
galaxies measured several times. This will be called the 
"auto-cross-identification" . Then, the cross-identification 
with LED A galaxies may start. 

6.1. Auto-cross-identification 

The "auto-cross-identification" is performed using a hier- 
archical method in which we merge step by step the closest 
objects. The definition of the distance of two objects i and 
j is the following: 



fc=l 



\Pik ~Pjk\ 
2<7 k 



(6) 



where, n p is the number of parameters (coordinates, di- 
ameters magnitudes...) for a given object, pik is the k- 
th parameter of object i, is the standard error of the 
k — th parameter. When two objects are merged they are 
replaced by a single one, whose parameters are the means 
of both. The final result does not depend on the order the 
original file is read. Note that, special care must be taken 
for periodic parameters, Right ascension and position an- 
gle (e.g., p. a. = Odeg is identical to p. a. — 180 deg; thus, 
the mean of p. a. = 3 deg and p. a. = 177 deg is p. a. = deg, 
not p. a. — 90 deg). 

The adopted exfc are given in the following table: 



Parameter 


Ok 


Ck2000 


30// 


^2000 


30// 


logD 


0.08 


logi? 


0.08 


It 


0.3 magnitude 


P 


5 / log R degree 



The standard error of the position angle is taken as a 
function of logR because its meaning vanishes for face-on 
galaxies. 

Objects are merged for dij < du m it, diimit beeing cho- 
sen from the distribution of all distances (Fig.||). By its 
definition, dij has the meaning of a Student's t-test dev- 
ided by 2 (it is thus dimensionless) . A la criterion corre- 
sponds roughly to du m a = 0.5. However, the value of a at- 
tached to each parameter is somewhat arbitrary, so is the 
definition of du m it- We adopted du m it = 0.55. This choice 
is guided by the minimum observed in the histogram of 
dij (Fig®. 

During "auto-cross-identification" a provisional DE- 
NIS number (called RED) is given to each entry of the 
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Fig. 6. Distribution of "distances" dij between two DENIS 
objects (i and j). This graph is used to chose du m it in the 
auto-cross-identification phase (see text). 
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Fig. 7. Distribution of "distances" dij between DENIS and 
LEDA galaxies (i and j). This graph is used to chose du mit in 
the cross- identification phase (see text). 



catalog. When a galaxy appears several times in the cat- 
alog, each original set of measurements is identified with 
the same RED number. Each entry will be cross- identified 
independently with LEDA galaxies, this will allow us to 
check the reliability of the cross-identification. 

6.2. Cross-identification with LEDA galaxies 

From the previous step we get an intermediate catalog in 
which each galaxy has a provisional number and all its 
astrophysical parameters (0:2000, ^2000, logDj, logRi, It 
and (3). Each entry of this catalog must now be cross- 
identified with LEDA galaxies in order to identify those 
already known. 

LEDA galaxies have similar parameters 02000 > ^2000, 
logD2s,, logB,25, Bt and 0, but diameters and magni- 
tudes are defined in the photometric B-band. The cross- 
identification is done by calculating the distance (in the 
mathematical sense, as defined by equation |^) between a 
DENIS and a LEDA galaxy. The coincidence is accepted 
if the distance is smaller than du m it- From a histogram 
of all distances between LEDA and DENIS measurements 
(Fig.^), we adopted di imi t = 1.0 which corresponds to the 
first minimum of dij ; -histogram (a pure Student's t-test 
would have given dumit = 1.29 at a 0.01 probability level). 
This limit is voluntarily conservative (i.e., small) because 
we prefer to miss a cross-identification than to merge two 
distinct galaxies. 

Because a given galaxy is cross-identified each time it 
appears in the catalog, 15945 were cross-identified sev- 
eral times with their original parameters. We reject 1881 
galaxies identified with more than one LEDA galaxy. The 
final catalog contains now 36247 galaxies. 



6.3. Last cleaning 

Many objects were kept in the catalog despite the fact 
that they were labelled "undefined objects" by the DA. 
They were kept because a known galaxy was very close. In 
the present release we removed these objects which nature 
may be questionable without further inspection. Indeed, 
some very faint galaxies in LEDA have only a few param- 
eters, so that the cross-identification is based mainly on 
coordinates. In some crowded field (clusters of galaxies or 
low galactic latitude) the accuracy of coordinates does not 
allow an identification secure enough. 10696 such objects 
were then rejected, leaving 25551 galaxies. 

Finally, we rejected galaxies with uncertain coordi- 
nates so that only coordinates based on the GSC reference 
are used. So, 5291 galaxies are rejected in the present ver- 
sion, leading to the final catalog of 20260 galaxies. These 
drastic rejections aim at maintaining a high quality level 
for this first catalog. In order to judge the quality more 
quantitatively we will now compare with other sources of 
data. 

7. Comparison with Mathewson et al. samples 

7.1. Magnitudes 

Magnitudes are calibrated by comparison with the mea- 
surements in I-band photometry made by Mathewson et 
al. (1992, 1996) which gives access to 2441 galaxies. This 
comparison allows us to correct for a possible variation 
of the zero-point. This variation has been explained by 
seasonal variation of the mean temperature of the cam- 
era El. Figure H shows such a variation described by the 
parameter C: 



the camera is now equiped with a regulated cooling system. 
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Night number 


C 


1850 - 2150 


+0.05 


2151 - 2350 


-0.08 


2351 - 2600 


+0.03 


2600 - 2800 


-0.05 



Note that the "night number" is a logical number, not 
a real night number. Each time the system is initialized 
the night number is incremented. This explains that after 
one year of running survey there are 2500 logical nights. 




I(Mathewson)- I(DENIS) 



2000 2200 2400 2600 

night number 



2800 



Fig. 8. Zero-point variation obtained by comparison with 
Mathewson's I-band photometry. 



An airmass correction is adopted using a typical value 
AJ/AsecC = 0.04+0.02. A check is made to control that 
there is no airmass residual. The residual is smaller than 
0.01 magnitude. The DENIS I-magnitude is then: 



I(DENIS) = 

71 

-2.5log(J2 fij - fbg) + 24.01 - 0.04secC + C (7) 

i=l 

The zero-point distribution of I (Mathewson) — 
I(DENIS) is Gaussian (Fig|]) with a standard deviation 
of 0.2 magnitude. If we assume that the error is identical 
for both systems the mean error on DENIS extragalactic 
I-band magnitudes would be about 0.14 magnitude. Be- 
cause the uncertainty on Mathewson et al. data is proba- 
bly smaller, the uncertainty on our I-band magnitudes is 
about ~ 0.18 magnitude. 

In FigM the comparison between I(Mathewson) and 
I(DENIS) is shown for galaxies with secure identification 
and being neither "multiple" nor "truncated" . The direct 
regression is: 



I (Math.) = 1.05 ± 0.02I(DENIS) - 0.54 ± 0.24 



(8) 



Fig. 9. Zero-point distribution after a tiny seasonal correction 

with the following standard deviation, correlation coeffi- 
cient and number of objects: a — 0.21, p — 0.977 + 0.004, 
n = 163 after 11 rejections at 3<r. 

Stricto sensu, the slope is not significantly different 
from 1, and the zero-point is not significantly different 
from 0. So, we are keeping the conservative solution: 
I(Math.) = I(DENIS). Among the 10 rejected galaxies, 
8 can be explained by localized poor photometric condi- 
tions (because they correspond to a loss of flux) . The mea- 
surements of corresponding nights will be counted with 
half weight. Two nights were rejected (night 2475 and 
2159) because they give rejections in the comparison of 
different photometric parameters. 

7.2. Diameters ans Axis Ratios 

Diameters and axis ratios are also compared with those of 
Mathewson ct al. samples. These comparisons are given in 
Fig.[ll] and FigJlJ, respectively. The direct regression are: 

logD(Math.) = 

0.96 ± OMlogD(DENIS) + 0.04 ± 0.04 (9) 

with a = 0.10, p = 0.88 ± 0.02, n = 170 after 4 rejections 
at 3(7. 

For the axis ratio, it is better to force the intercept to 
be zero as generally admitted (see de Vaucouleurs et al., 
1976). This avoids to have negative axis ratio after the 
application of the regression. The result is thus: 

logR(Math.) = 0.94 ± 0MlogR(DENIS) (10) 

with a = 0.10, p = 0.85 ± 0.02, n = 172 after 2 rejections 
at 3(7. 

None of these regressions is significantly different from 
the absolute identity. So, we will keep: logD(Math.) = 
logD(DENIS) and logR(Math.) = logR(DENIS). The 
standard deviations are 0.10 for both logD and logR. 
Again, if we assume the same error on both systems the 
error on logD and logR is about 0.07. 
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Fig. 10. Comparison of extragalactic I-band magnitudes from 
Mathewson et al. and from DENIS 




logD(DENIS) 

Fig. 11. Comparison of I-band isophotal diameters from Math- 
ewson et al. and from DENIS. 




0.5 



logR(DENIS) 

Fig. 12. Comparison of I-band axis ratios from Mathewson et 
al. and from DENIS 

8. Comparison with LEDA galaxies 

8.1. Equatorial coordinates 

The coordinates are compared with coordinates given in 
LEDA. Only coordinates known as 'accurate' in LEDA 
(i.e. with a standard deviation less than 10 arcseconds) 
are used. The plot of AS and Aa.cosS in Figure [l^ 
shows that there is no systematic distorsion (A means 
(DENIS) - (LEDA)). The standard deviation is 6.5 arc- 
seconds and 6.7 arcseconds for a.cosS and 5, respectively. 
Assuming an error of the same amplitude in LEDA and 
DENIS coordinates gives an uncertainty of 6.6/v / 2 arc- 
seconds for DENIS right ascension and declination, i.e. an 
uncertainty of 6.6 arcseconds for the position of a galaxy. 

In fact, the plot of AacosS vs. Ai5 exhibited two ab- 
normal zones with a systematic shift of about 30 or 40 
arcseconds. This problem appeared near the zenith. Thus, 
only objects with coordinates based on the GSC stars are 
kept in the present version. In the final catalog, the coor- 
dinates will be obtained through a full astrometric solu- 
tion (mosaicing frames along each strip and with adjacent 
ones) and this problem will be solved without rejecting 
galaxies. 

8.2. Position angle 

Position angle is important for studies on the orientation 
of galaxies, but also for identification. However, for nearly 
face on galaxies, it becomes very uncertain. The compar- 
ison is made only for galaxies with logR(DENIS) > 0.5. 
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Fig. 13. Comparison of coordinates from LEDA and DENIS. 



Fig. 

NIS. 



14. Comparison of position angles from LEDA and DE- 



The result is given in (Fig.nih. The direct regression gives: 



P(LEDA) = 1.009 ± 0.003[3{DENIS) + 0.07 ± 0.26 (11) 

with cr = 2.65, p = 0.9987± 0.0001, n = 408 after 21 rejec- 
tions at 3cr. Among the 21 rejections, only 2 correspond 
to real discrepancies. 

The uncertainty on the measurement of (3 is 2.7deg.. 
This excellent agreement of position angles validates the 
reliability of our cross-identification. 



8.3. Photometric morphological type code 

It is interesting to obtain at least a very crude estimate 
of the morphological type code. This is particularly im- 
portant when we plan to start a HI follow-up for which it 
is compulsory to avoid elliptical galaxies. The log of the 
standard deviation of pixel intensities is significantly cor- 
related with the morphological type code extracted from 
LEDA. The relation is almost linear. However, it appears 
that the solution depends on the size of the galaxy (all 
very small objects appear identical). The slope and the 
intercept are found linearly correlated with the log of the 
number of pixels. In the scale between -5 and 10 defined 
in the Second Reference Catalog of Bright Galaxies, the 
photometric morphological type code is calculated as: 



Where o{fij) is the standard deviation of flux of all pixels 
of the matrix representing the galaxy : 



with 

A = 



(n-l) 



1 



-0.0187 logn- 0.0183 



B 



1 



0.2774 log n+ 1.7846 



(13) 



(14) 



(15) 



T(DENIS) = A logK/y)) + B 



(12) 



where n is the number of pixels. 

The comparison of codes T(DENIS) and T(LEDA) 
is given in Figure |l5|. The correlation is clearly significant 
(the correlation coefficient is p = 0.69 ± 0.01) but the 
standard deviation is large (a(T{DENIS)) = 2.5). This 
allows us to classify galaxies into 'Early', 'Intermediate' 
and 'Late' types, with no finer subdivision. In the cata- 
log, photometric morphological type smaller than —5 or 
greater than 10 will be set to —5 or 10, respectively. 

9. The catalog 

9.1. Description of the catalog 

The catalog of galaxies which is described here is the first 
of a series which will be released during the progression of 
the survey. At the end of the survey a deeper and complete 
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Fig. 15. Comparison of morphological type codes from LEDA 
and from the I-band photometry. 

catalog will be produced at PDAC. The present one, is 
necessarily limited to the first observations (one year). 

The catalog contains 20260 galaxies. Among them 
14518 are new galaxies while 5742 are galaxies already 
present in LEDA. 

Each galaxy is numbered with a provisional inter- 
nal number (RED, for Rapid Extraction from DENIS), 
.and an identifier from the LEDA database (PGC/LEDA) 
□. Galaxies are also identified with an alternate name 
taken in the following catalogues: NGC, IC (Dreyer, 
1988), UGC (Nilson, 1973), ESO (Lauberts, 1982), MCG 
(Vorontsov-Velyaminov et al., 1962-1974), CGCG (Zwicky 
et al., 1961-1968), IIZW, IIIZW, VIII from the cata- 
log of compact and eruptive galaxies (Zwicky, 1971), 
IISZ (extension of the same list by Rodgers et al., 1978, 
FAIR (Fairall, 1977-1988), HICK (Hickson, 1989), KUG 
(Takase and Miyauchi-Isobe, 1984-1989), IRAS (IRAS), 
MK (Markarian, 1971-1977 ), UM (University of Michi- 
gan list, MacAlpine et al., 1981). The detailed references 
to these catalogs are given in Paturel et al. (1989). Two ad- 
ditional catalogs are included: Saito et al. (1990-1991) and 
Dressier (1980). The notations for corresponding galaxy 
names are SAIT (e.g. SAIT 69-1, for object 69 in list 1 of 
Saito et al.) and DRCG (e.g., DRCG 39-41, for galaxy 41 
in cluster 39, the numbering of clusters is made according 
to the order in the original publication). 

For each galaxy, the catalog gives the weighted mean 
parameters obtained as described in the previous section. 
Actual mean errors a are calculated according to Paturel 
et al. (1996) using individual mean errors deduced in 
previous section. For I-band magnitudes the mean error 
is taken as a function of the magnitude itself. An esti- 

2 connection: 

http:/ /www-obs. univ-lyonl.fr/leda 
or 

telnet: lmc.univ-lyonl.fr logindeda 



mate made from Fig. [H] gives:cr(I) = 0.05B T - 0.51. The 
weights used for calculating the weighted means are the 
inverse square of actual mean errors. For nights suspected 
to be of poor_photometric quality the actual mean error is 
divided by 30. Further, when the quality of the matrix is 
uncertain ('peculiar', 'multiple' or 'truncated'), the weight 
is divided by 2. This correcting coefficient is deduced from 
the comparison of I magnitudes with poor quality objects. 

The quality of each individual matrix is coded as fol- 
lows: 'Normal' galaxies are coded as 1, 'Peculiar' galaxies 
as 10, 'Multiple' as 100 and 'Truncated' as 1000. This code 
is added for each measurement. So, the resulting code gives 
immediately the number of independent measurements 
(sum of digits) and the quality of each of them. For in- 
stance, the code 1002 means that there are one truncated 
matrix and two normal ones (i.e., 3 independent measure- 
ments). 'Truncated' means an overestimated magnitude, 
on the contrary, 'Multiple' means that the magnitude of 
the considered galaxy is underestimated. 

The different columns are the following: 

— Column l:Provisional DENIS identifier. This identi- 
fier is just an internal number given for easy identifica- 
tion. This number will be replaced by an appropriate 
final DENIS number. 

— Column 2:LEDA identifier. This number allows rapid 
retrieval in the LEDA database. This number is iden- 
tical to the PGC number for LEDA number less than 
73197. 

— Column 3:Alternate name according to a hierarchy: 
NGC, IC, ESO, UGC, MCG (see text above). 

— Column 4:Equatorial coordinates for epoch 2000, in 
hours, minutes, seconds and tenths for the Right As- 
cension and in degrees, arcminutes and arcseconds for 
the declination. The actual mean error on the position 
of the galaxy is given in the same column, next line in 
arcseconds. 

— Column 5: Total I-band magnitude I(DENIS) with 
its mean error on the second line. 

— Column 6: Isophotal diameter in log scale 
logD(DENIS), where D is in 0.1 arcminute follow- 
ing the notation of de Vaucouleurs et al. (RC1, RC2). 
The actual mean error is given on the second line. 

— Column 7: Isophotal axis ratio in log scale 
logR(DENIS), where R is the ratio of the major to 
the minor axis. The actual mean error is given on the 
second line. 

— Column 8: Position angle (3{DENIS) in degrees, 
measured from North Eastwards. Its value lays 
between and 180 degrees (excluding 180). For 
logR(DENIS) < 0.5 the meaning is poor. This is re- 
flected by the actual mean error given on the second 
line. 

3 These nights were detected by a 3<r rejection rule, so one 
can admit that their mean error is three times the typical mean 
error. 
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— Column 9: Photometric morphological type code 
T (DEN IS). The coding is given according to RC2 
(i.e. T from —5 to 10 for Elliptical to Irregular galax- 
ies). The mean error is given on the second line. 

— Column 10: Quality code of the object (1000=trun- 
cated object, 100=multiple object, 10=peculiar object, 
l=normal galaxy). The sum of digits gives the number 
of independent measurements. 

A sample page is given at the end of the paper. The 
catalog is available in electronic form at PDAC and Lyons 
Observatory. 

9.2. Distribution on the sky 

The distribution on the sky is given in equatorial coordi- 
nates for the epoch 2000 (Figure |l6|). The strips are clearly 
visible. The zone called 'Equatorial' (5 between +2 deg 
and —28 deg) is not so well covered because observations 
are avoided in this zone when Moon is bright or when the 
wind is strong. Further, many frames are rejected in this 
zone at the end of the strip (near —28 deg) because of the 
shift in header-coordinates. No attempt is made to reach 
low galactic latitudes. This is done independently in J and 
K bands. 

9.3. Completeness limit 

If we assume a homogeneous distribution of galaxies (N cx 
r 3 ), the plot of the number of galaxies (in log scale, i.e. 
logN) brighter than a given magnitude limit Iu m it versus 
this magnitude limit should be linear with a slope of 0.6. 
This plot (Figure |T^ ) is quite linear up to Iumit = 14.5 
with a slope of 0.62 ± 0.02. The sense of this completeness 
limit must be explained. It means that in surveyed direc- 
tions all galaxies brighter than 14.5 are detected. Because 
the sampling is made randomly and 1/4 of the survey is 
presented here, the completeness limit in apparent mag- 
nitude of this first DENIS catalog is Iumit = 14.5. 
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Fig. 16. Distribution of the 20260 galaxies of the present DENIS catalog presented with a Flamsteed's equal area projection 
in equatorial (2000) coordinates. 
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TABLE 1: Sample page of the present catalog. 



DENIS-RED 


LEDA/PGC 


Alternate Name 


R. A. DEC. (2000) 


I 


logD 


logR 


p. a 




T 


code 


(1) 


(2) 




(3) 


(4) 


(5) 
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0. 
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95 


9 
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4 


.21 




.052 
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21 


5 
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6 


.36 




.099 
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6 
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